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BOUNDARY VALUE AND EXPANSION PROBLEMS OF ORDINARY 
LINEAR DIFFERENTIAL EQUATIONS* 

BY 

GEOKGE D. BIEKHOFF 
Introduction. 

Let p 2 (x), p 3 {x), •••, p n (x) be functions of the real variable x on the 
closed interval ( a , b), which are continuous with their derivatives of all orders. 
Write 

T . s d n z , . d n ~ 2 z , s 

(1) 

M(z)^(-iy~+ * +(-iy- t ^=f[p i (x)z]+.-.+p n (x)z. 

With the linear differential equation of the nth order in u 

(2) _Z/(M) + \w = 0f 

and n linear homogeneous conditions in u(a), w'(a), •••, w (n_1) (a), w(6), 
«'(&), ■--,«<— ^(ft), 

(3) FiOO-O, W 2 (u) = 0, -.., W n (u)=0 
we associate the like adjoint differential equation 

(4) M(v) + \v = 
and n like adjoint conditions 

(5) ^(^ = 0, V 2 (v) = 0, -.., V n (v) = 0. 

For certain characteristic values of the complex parameter \ there will exist 

*The second part of a paper presented to the Society (Chicago), March 30, 1907, under a 
different title. The first part of this paper has been printed on pages 219-231 of this volume. 
Received for publication May 12, 1908. 

t It is not an essential generalization to write instead of (2) 

a(») £ + ft(*)^! + •■• + *.(*)« + ^(*)»=0 

if j) + 0, P+0, and p /P is real. We therefore restrict (2) as stated. Max Mason has 
treated a special case n = 2 with the restriction .P+0 removed, by different methods : these 
Transactions, vol. 8 (1907), p. 427. 

Trans. Am. Math. Soc. 85 373 
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a solution* u(x) of (2), (3) or v(x) of (4), (5). These values are the same in 
both cases : let them be 

\ ' ^2 ' 

with corresponding solutions 

M i(aO> u 2 ( x )> •■■> 

v^x), v 2 (x), •••. 
These functions u t (x), v.(x) are such that 

(6) I u i (x)v J (x)dx = (*+i)- 

This property leads us to the formal expansion of a given function J\x) on (a, 6), 

» I y(»)« 1 (*)*» 

How are the characteristic values distributed in the X-plane? What is the 
nature of the solutions u { (^x) , v t (x)? In what sense does the expansion repre- 
sent f(x)? These are the questions considered in this paper. 

We begin with the derivation of the formal properties of the boundary value 
problem (§1), and of the expression for the sum of n terms of (7) by means of a 
contour integral (§ 2). 

There follows the more intimate study which is based on certain facts concerning 
the asymptotic nature of the solutions of (2) and of (4) when | X | is large. These 
facts are derived as an application of my paper in a preceding number of these 
Transactions (§3). The distribution of the numbers V and the nature of 
the expansion is then obtained (§4 and § 5). Finally the contour integral is 
evaluated and the representation theorem proved (§ 6). The expansion is found 
to behave like a Fourier series except in the vicinity of x = a and x = b . 

Liouville was the first to introduce the notion of adjoint conditions in a special 
case f and to consider the related expansion. The results of the present paper are 
known for the real self-adjoint case n = 2 . % Westfall has proved a repre- 
sentation theorem for the real self -adjoint case n = 2m, providing that /"(as) 
and its first n derivatives are continuous. § 

* By a solution we mean always a solution not identically zero. 

■fLiouville's Journal, ser. 1, vol. 3 (1838), p. 561. Professor E. H. Mooek suggested 
to me the possibility of generalization. 

X A. Kneser, Mathematische Annalen, vol. 58 (1901), p. 81 ; Dixon, Proceedings 
of the London Mathematical Society, aer. 2, vol. 3 (1905), p. 83. 

2 W. D. A. Westfall, Zur Theorie der Integralgleichungen, A. Myllek, Qewohnliche Diffe- 
rentialgleiehungen hoherer Ordnung in ihrer Beziehung zu den Integralgleichungen. Gottingen 
dissertations (1905 and 1906). Reference should also be made to Hilbebt's antecedent papers on 
integral equations in the Gottinger Nachriohten for 1904. 
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§ 1. Th,e formal nature of the boundary value problem. 

We must first make precise the adjoint conditions referred to. 

Definition. Let W l (u), • ••, W K (u) be n given linear forms in u(a), 
«'(«), ...,«f-"(o),«(6), u'(6), ■■■,u^\b) and W n+l (u),..-, W 2n (u) 
be any n further linear forms so chosen that TF, , • • • , W 2n are linearly inde- 
pendent. Then in the identity 

(8) CzL(y)dx = P(y,z) ""* + C yM(z)dx* 

where -P(y, s)|^r* is a bilinear form in y(a), y'(a), ■••, y (n_1) (a), y(b), 
y'(b), ...,y*-»(b)zndz(a),z'(a),...,z^(a),z(b),z'(b),...,z^(b), 
we can write 



(9) P{y,*) 



= E^(y)r*_<(*). 



in which V t (z) are linear in z(a), z'(a), ■ ■ ■ , z ( "- 1) (a) , 2(6) , 2'(6) , • • • , z (n -^(6) 
and linearly independent. Then Fj(v)=0 (£ = 1,2, •••,?!) are the adjoint 
conditions to the given conditions W.(u) = (£ = 1, 2, • • •, n). Any set of 
conditions W { (u) = (£ = 1, 2, •••,«) equivalent to W t (u) = (£ = 1, 2, • • •, n) 
by linear combination, with any choice of W n+l (u), ■ •-, W 2n (u) leads to a set 
of adjoint conditions V t (v) = (i=l, 2, ••-, n) equivalent to V t {v) = 
(i = 1, 2, ••.,«) by linear combination. 

Conversely, given ^(-y) = 0(i = 1, 2, •••, w) we choose V n+l , • • •, P^ n as 
above and find JF i (i«) = 0(£ = l, 2,--«, w) to be the adjoint conditions. 
Hence if the problem (4), (5) is adjoint to (2), (3), so also is (2), (3) adjoint to 

(4), (5). 

The properties which this section proves are stated in I— III. 

I. If for \ = \* a solution u*(x) of (2), (3) exists, a solution v*(x) o/"(4), 
(5) will also exist for X=\* ; ifu*(x) is unique (except for a constant factor), 
v*(x) is also unique (except for a constant factor). 

Let u* be the given solution of (2), (3) for \ = \*. Then we have 

W l (u)= W 2 (u*) = ...= W n (u) = 0, 
and for some j 

as otherwise we should infer 

u'(a) = u*'(a) ■. u^\a) = u\b) = u\b) = . • • = u^ l \b) = 0. 

In the ?i-fold linear spread of solutions v(x) of (4) at \ = \* there will be at 
least one, say v*(x), which satisfies the n — 1 linear homogeneous conditions 

V n _ i (v) = {i + j). 

* CI. Schlesingkb : Handbuch der Theorie der linearen Differentialgleichungen, vol. 1, p. 54, 
formula (I). 
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Substitute in (8) y = «*, z = v* and we obtain 

-X* f vudx = W n+J (u)V n _j(v m )-\' f uv'dx. 

Therefore V n _.{v) = also, and v* will satisfy (4), (5). 

If u is unique, v* is also unique. To prove this assume if possible that 
there were two linearly independent solutions v* and «** satisfying (4), (5) when 
X = X*, while u* is unique. Then we could choose i and j so that 



4=0; 



V n+i (v*) V n+J (v) 

V n+i (v<*) V n+J (v") 

otherwise we should have constants c, d, not both zero, such that 

V n+l (cv* + dv**) = (i = l, 9, •■•,«), 

and v = cv* + dv* would fulfill the 2n conditions 

V,(v) = V 2 (v)=... = V 2n (v) = 0, 

which is not possible since the V 1 , •••, V 2n are linearly independent. Choose 
now u* linearly independent of u to fulfill the n — 2 conditions 

W^u**) = 0, (Z 4= ji — i or n —j) 

where I = 1 , 2 , ■ • -, n. Writing in (8) y = u**, z = v* and y = «**, z = a**, we 
obtain, on simplifying, 

T^(«") F n+< (**) + TT„_,( «**) F„ + .(^) = , 

TF^(0 F^(0 + TF^(«-) F„ + ,(0 = . 

From this we deduce that also 

Therefore u m would satisfy (2), (3). This is impossible. Hence if u* is unique, 
v* is also unique, which we were to show. 

Definition. If X = X* is a characteristic value of X for which one and but 
one solution of (2), (3) and (4), (5) exists, X* is said to be a simple characteristic 
value. 

II. If ' y x , y 2 , • • ■ , y n are n linearly independent solutions of (2) at X = X", 
the condition that X* is a characteristic value is that the determinant 



A = 



^(3/,) W 2 {y 2 ) ... W 2 (y n ) 



wm W n (y 2 ) 



WM 
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vanishes; the condition that A* is a simple characteristic value is that not all 
the first minors vanish. 

The general solution of (2) is 

u = c l y l + c 2 y 2 +--.+c n y n . 

Under the condition A = at least one determination of c x , c 2 , • ■ • , c n is possible 
such that 

TT 1 ( M ) = 0, W 2 (u)=Q, ..-, W n (u) = 0. 

This determination is unique when not all the first minors of A vanish. 

III. Jfu t (x} and v.(x) belong to distinct characteristic numbers \. and A.., 
then is 



£ 



u^.dx = 0. 



For write y = u t , z = v. in (8). Because of (9) we deduce 

X6 fib 

v j u i dx= — \. I u { Vjdx, 

whence the desired relation follows. 

In important particular cases the problem (4), (5) may be precisely (2), (3) 
with v in place of u. We then say that the problem (2), (3) is self-adjoint. 
Or (4), (5) may be precisely (2), (3) with v in place of u, except that M (v) is 
the negative of L(v ) . We then say that the problem (2), (3) is anti-self-adjoint. 

§ 2. Expansion as contour integral. 

It is easy to prove that when \ is not a characteristic number there exists a 
unique G(x, s; A.)* such that the solution <£ of 

Z(4»)+X<£ = a,, W^)= W 2 (<f>)=---= W n (4>) = 0, 

is given by 

<£= 1 G(x, s; \) (o(s)ds; 

and likewise that there exists a unique JT(x, s; X) such that the solution <f> of 

M(4>) + \4> = *, r l {<t>)= v 2 (<t>) v n (<t>) = o, 

is given by 

<f> = I IT(x,s; \)m(s)ds. 

* Professor Bocher has defined the Green's function Q for an ordinary linear differential 
equation of order n tinder a special form of conditions, and has stated the principal properties in 
the Bulletin of the American Mathematical Society, vol. 7 (1901), p. 297. 
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The explicit formula for G is 



(10) G(x,s; X) = 



N(x, s; X) 
~"A(XJ-~ 



y t (x) y 2 (x) ••• y n (x) G(x,s;\) 

WM W 2 (y 2 ) ■■■ W 2 (y n ) W 2 {G) 
WM W H (y 2 )-- W n {y n ) W n (G) 



(-1)" 



in which 






(ii) 



£(«, s; \)='±i- 



^(yj W n {y 2 ) ... W n (y n ) 

Vi( x ) y 2 ( x ) ••• y„0) 
#-«(«) 2/ 2 »- 2, ( s ) ••• yT 2) ( s ) 



yi( s ) y 2 ( s ) 



y»(«) 




yi( s ) y 2 («) ••• y„(«) 

(+ if a;> s, — if a; <s). 

Here y 1 , y 2 , • • •> y„ are any n linearly independent solutions of (2). There is a 
like formula for ff(x, s; \). Furthermore 

G(x,s; X)=H(s,x; \).* 

The function G is analytic in X (since y t , y 2 , • • • , y n may be taken analytic 
in \) except for a possible pole when A (A.) = 0, i. e. by II, § 1, wAew X is a 
characteristic value. 

If\ = \.is a simple characteristic number for which G has a pole of the 
first order, the residue is 

u i( x ) v i( s ) 



where 

Write 

(12) 



j u.(x)v.(x)dx 

I u i (x)v i (x)dx =J= . 

H(x, s) 
G(x, s; \) = x _\ + <r(cc, s; \), 



* For formulas like these and their proofs see Westf all's dissertation, \\ 6, 7. 
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where <r(x, s ; X) is analytic at \ = \. In view of (10) we have 

JV(x ,s; X.) 

s ( x ' 8 ) — a 7 ^ — ' 

where JV"is continuous with its first n derivatives since the coefficient of G is zero. 
Further N satisfies (2), (3) for all \ so that M(x, s) must be a solution of 
(2), (3) in x. Because of the relation between G and H, we infer also that 
H(x, s) is a solution of (4), (5) in s. By the definition of a simple character- 
istic value it follows that the residue is 

M(x,s) = c i u.(x)v i (s).* 

It remains to determine c r 
By (12) we have 

lim [(X — \) G(x, s; X) — c.u i (x)v i (s)'] = 0, 

A=A f 

whence 

lim|(X — \) I G(x, s; \)u i (s)ds — c^^x) I v.(s)u.(s)ds =0. 

A=Aj L *J a *Ja _| 

But since 

£( Wj ) + \ w< =(\-VK, W l (u i )=W 2 (u i ) = ...= W n (u i ) = 0, 
we have by the fundamental property of G 

(\ — V) I G(x, s; \)u.(s)ds = u.(x). 

%J a 

Substituting this value above we find 

lim ^(cc) — c t u.(x) I v i (s)u.(s)ds = 0. 

A=A| L Ja J 

i J^v i {s)u i (s)ds=l, 



Therefore 

-.6 

c.- 



out of which c. is determined ; this proves that the residue has the stated form. 
If r be a contour in the X-plane which encloses \, X 2 , • • •, X n , we conclude by 
the above that 



\ r r" *=, I f( x ) v i( x ) dx 
7== #(*,*; X)/(*)<Mx = £^- „ 4 ( 



* In this proof certain points of logic are obviously slurred over. 
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This is the sum of n terms of the formal expansion, providing that \, \ 2 , • • •, X n 
are simple characteristic numbers at which G has a pole of the first order. At 
values \. for which these conditions are not satisfied, the corresponding term 
of the formal expansion (7) is to be replaced by 



j l ' B.(x, s)f(s)ds 



where H t (x, s) is the residue. Thus in all cases we have an expansion each 
term of which is uniquely determined. We omit the development of H { (x, s) 
in these more complicated cases. 

§ 3. The solutions of L(u) -f- Xu = and M(v) + Xv = when \\\ is large. 

As a preliminary to the deeper study of the problems which have been pre- 
sented, I apply the results of my paper printed on pages 219—231 of the present 
volume to the differential equations L(u) + \w = and M(v) + Xv = 0. 

If we write X = p", the first of these differential equations can be written in 
the form there treated [see (1), loc. cit.] , 

d n u , .d n ~ x u . , 

^ + /»_,(», p)^=v + ■•■ + PXK P) u = °» 

a n _ l{x ,p) = 0, a n _ 2 (x,p) = ^, -.., a (x,p) = l + ?j^. 

The coefficients obviously satisfy the restriction imposed, and the equation for 
w t (x) [see (4), loc. cit.], is 

(14) < + l = 0. 

Thus wjj, w 2 , ■ • •, w n are constants. We also find (proof of Lemma 1, loc. cit.) 

(15) u m (x) = l (< = i,8, ••,»)• 

In order to state in explicit form the theorem for this case, it remains to con- 
sider the regions S (definition, p. 220). Inasmuch &s w x ,w 2 , • • •, w n are con- 
stants we can for each p choose the indices 1 to n so that 

j?(p«0 = -»0«°,) = ■■■ = -» o»o 

for every x on (a, b). Hence every p is on some region S, = argp = yjr. 
At the bounding rays the ordering changes so that for some i and j 
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R(pw i ) = B{ P w.) 
An easy computation shows that the regions S are the regions 
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(16) 



^argp^^ 



(*=(), 1, .... 2n-l). 



n ° ' n 

The differential equation M(v) + \v = also is of the same type when we 
put \ = p n . The roots w t (x) are 

(17) — w x , — w 2 , •••, — w n 
and in the formal solutions we have 

(18) u. (x) = l (t = l, 2, •••,»). 

The regions S are the same as for L(u) + \w = 0. 

Let i? be a generic notation for functions of p (and other variables), bounded 
when j p | is large. The application of the theorem referred to gives : 
On any region of the p-plane 



(19) 



# : — = arg p ~ — 



(20) 



there exist n independent solutions, 

2/i> y 2 » ■••» yn / x ( M ) + p" M = ' 

analytic in p and such that on this region 

y i =u i {x,p) + e^-^- p l, 

dy- d „ , , E. 

Ix=dx U ^^ + ePW ' {X ~ a) ^ 



(21) 



d n ~ x y. d"- 1 , s , s E. 



dx n 



(22) 



where 



dx' 

z . = v .( x ,p) + e -^-^—\ 

dz. d , „ , ^ E, 



dx dx ' 



1 z Jn-l j^ 

n-1 — ,7,v,n-l V i\ X > ") T e „m-n+H 



<fe 



c?x n 
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«,(», p) = *■«««> [l + ^ + • • • + ^^1 , 
(23) 

ei ( a> ,p) = e-P-K—)[l + ^ + ... + ^] 1 

Wj , w 2 , • • • , w n being the n roots o/" io" + 1 = , m any positive integer. 

§ 4. Distribution of the characteristic values. 

The condition that A* is a characteristic number is that the determinant A 
vanishes (II, § 1). The i-th. element in the j-th column of A is W t (yj) where 
y x ,y 2 , • • •, y n are any n linearly independent solutions of _£(«) + A* w = 0. 
In order to treat the equation A = we take y t to be the y i of (20) § 3 where 
\ = p". In addition we assume the conditions W- to be normalized as follows. 
Keduce the number of conditions W-(u) = of order n — 1 [i. e., containing 
either u (n ~ 1} (a) or m ( " _1) (6)] to a minimum, at most 2, by linear combination. 
Then, in those that remain reduce the number of order n — 2 to a minimum, 
at most 2, by linear combination. Continue in this way as long as conditions 
remain. The normalized conditions will have the form 

(24) W i (u) = W ia (u) + W ii (u) = 0, 

W ia (u) B a jtt <«(a) + g a..u^(a), 

ir («)-/s 4 «w(&) + £\«y>(&), 

in which no three successive &'s are equal. 

Definition. Let the w. be taken in every order such that for some p 4= 

R(pw l )<B{pw 2 )<.--<R{pw n ). 

If n = 2fi — 1 and always neither $ = nor # x = , where O and X are defined 
by the identity 

a x w\ l •■■ «! w*«_i (ai + sA)w* 1 frwJiVi "• ^i M » 



(25) *„ + '.* = 



«»m{" ••• a„< n -i (a„ + s/3„X» &m*Vi •" &> w *" 



or if ra = 2/u. and always neither ^ = nor 2 = , where O , lt and 2 are 
defined by the identity 
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(26) d + 9 lS + -f = 

£-1 (« 2 +«ftW («2+7^)<+l &«£+ 



+2 •' 



a 2 ixfy • • • a 2 wZ 



&«£ 



&«£ 



.x (a»+«&)«£ (« B +jfl.)«£n &«£+« ••• A 



<*„wT 



•V^"- 






the conditions TF)(m) = (i = 1 , 2, •••, w) are said to be regular.* 

We now proceed to prove the 

Theorem. Jf the conditions W i (u) = ( i = 1 , 2 , • • • , w ) are regular, the 
characteristic values of \ for (2), (3) and (4), (5) are in general simple f ana" 
ybnw a pair of infinite sequences \ n , \ n , (1 = 1, 2, • • •) smcA £Aa£ 



(27) 






where g v , g Uj are constants and | E lt | < Jf, | i? IK | < M. 

We consider the equation A = in the p-plaue, fixing attention on some region 
(19) for I = l . The transformation \ = p n makes two adjacent regions S 
correspond to the entire \-plane. It is convenient to consider separately the 
cases n = 2/x. — 1 and n = 2/x. . The starting point is the set of relations 
[deduced from (21)] 



(28) 



[m—\ - Tp "I 

[m-1 Q W ~\ 

A + g^f'+^rJ. 



* This apparently gives one condition for each region 8, but these reduoe to two for n = 2/i — 1 
and to one when n = 2ft . It is worthy of note that the Stubm-Liouville boundary conditions 
( n = 2 ) of the form 

hu{a)+ku'(a)=0, lu(b)+mu'(b)=0, 

and also the periodio boundary conditions ( n = 2 ) 

«(a)=u(&), u'(a)=«'(6), 

are regular. An example of non-regular conditions ( to =2 ) is 

u(a)=0, « / (a)=fc«(i). 
tSee§l. 



384 G. d. birkhoff: expansion problems [October 

Case 1. n= 2/* — 1. 

By definition of the regions S we can choose indices 1 to n so that 

Ripw,)^ R{ P w 2 )^ ■ ■ .S H( P w n ) 

for the region (19) under consideration. Making this choice we find readily that 
on S 

B( P w i )<0 (i = i, 2, ■■;it-l), 

(29) 

B(pw.)>0 (i = f + l, ^ + 2, ■■-,»), 

while 

(30) ^(pwJ-0, argp = ^L±I),r. 
From (28) and (29) we see that in the formula 

^(%) = W^y,) + W a ( yj ) = 1, 2, ...,^-i), 

the term W ib (y J ) can be absorbed into the term W^y.) since e^ 6- ^ is small 
of an infinite order in 1/p on #; that is, we can write 

(31) W i (y j .) = ( P w i .)^a i + X,^ + -^ W =l l8| ... ,„-!). 

By like reasoning we find 

r m-l o 77 n 

(32) TF^y,) = e^-^,)*!^ + g ^ + - J (i=^ + l, ^ + 2, ••-, »). 
Also from (28) 

[m— 1 _ E* "I 

-,+ g-^ + ^J 

r m-l O ITr "1 

If we substitute in A these values of W^y.) as given in (31), (32), (33), and 
remove the factors p ki from the i-th row (i = l, 2, •••, w), and the factors 
e pw,(b-a) f rom t ne j.th column (._; = /* + 1, /* + 2, • • •, »), i. e., altogether the 
non-zero factor 

n n 

n /' • n e n{b - a \ 

the condition takes the form 

(34) A=^ + E^+^J + e^->^ 1+ g^ + ^) = 0, 

n n _ 

A = IJ p** • J J gPWAft-o) . A . 
i=l i=n + l 
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The quantities and 6^ are the & and Q x of (25). In fact (25) was formed to 
give these terms. 

Since A was analytic in p, A will be also analytic in p. 

The expressions of the first minors of the /it-th column of A , 

b E 

cannot have all their first terms a = since then 6 = 6 X = , contradicting 
the hypothesis. Hence for large values of | p | not all these first minors vanish. 
If we return to A, we see that accordingly not all its first minors vanish. The 
characteristic numbers when | X | is large must therefore be simple. This was 
part of the theorem. 

From (34) we obtain the equivalent condition 

*"**-* = 8(p), B(p)= -% + '£% + £ 

u i i=\ P H 
by solving for e p, V (6 ~" a> . From this we infer 

(35) pw ll (b-a) = log8(p), 

log S(p) = log(-^+£^ + |± 2^1/^1, 

the final form of the condition, in which of course E is analytic in p on S. 

In view of this condition, it is obvious that the values of p on S which satisfy 
A = have when | p | is large the form 

w*V~-i log \ IJ ^h E 

where I is a large positive integer. Since by (30) 

arg « M =, _ ^_^ v± _ , 

we see from (36) that the =fc sign is to be so chosen that the p t approach asymp- 
totically (in an angular sense) the bisecting ray of S. Let now p describe a 
circle of fixed small radius r about 



2W^t lo g(-<r°) 
... / j. _\ + ... / j. ~\ • 



w^b—a) w^{b — a)' 
This will be wholly within S when | p | is large ; and 

arg [pm^b - a) — log 8(p)~\ 
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increases by 2w if the inequality 

\r(b-a)\> 



Up- ' r 



obtains as is true for | p | large. This proves the unique existence of />, when 
| p | is large, and the distribution of values p is determined. 
It remains to return to the X-plane ; /S becomes S' 

S': Z 7rSarg\^(Z + l)u. 

The values of X approach the positive or negative axis of imaginaries according 
as l is even or odd, and by (36) have the form stated in the theorem where for 
some , 6 l 

.iog(4;) 



9u or g ui = 



torV-1 



This completes the proof when n = 2fi — 1 . 

Case II. n = 2/x . 

The proof is analogous to that for Case 1 . If pw x , pw 2 , • • ■ , pw n have 
increasing real parts on S we find 



while 



B(pw i )>0 (i=^ + 2, ft + S, ---,»), 

JR(pw li )= H(pw IJi+l ) = 0, argp = either-^ or ° ^ " ", 



and v> +1 = — w M . The region S is of course defined by (19) as before. We 
now find W.{ffj) to be of the form (31) for J = 1 , 2 , • • •, fi — 1 and of the form 
(32) for j = j»+2,(t+3 ) '",»; also for j = fi, n + 1 we find 

1=1 P P J L *=i /> P J 

r m ~ x n PI 

w^ +1 )= (- /*v P [*<+ § ^ + £ J 

r m-X Q J? -I 

writing w> M+1 = — u^. 

These expressions W^y.) are substituted in A and a factor 

n n 

JJ.P ki ' II e^* - ^ 
is removed. The equation A = becomes then of the form 



1908] OF OKDINABY DIFFERENTIAL EQUATIONS 387 

/ m— 1 f, Tf \ / to— 1 J J? \ 

(38) A = (. + g^ + -) + e^-(. 1 + g^ + ^) 

(TO— 1 „ J? \ 



A=H/ 3 *' II e" w ' <b - a) -A, 

#„,#,, 2 being defined as in (26). The equation (38) is quadratic in e" w ^'- a \ 
This equation (38) will yield the forms stated in the theorem for the charac- 
teristic values \ u , \ m on S. They must lie near to one of the bounding rays, 
say arg p = (Z fl + 1 )tt/ti. If we consider the necessary form of the expansion 
on the other region /S adjoining this ray, i. e., on the region 

(*0+l)7T < < q + 2)7T 

— — = are o s ^-^ — 

n &r n 

we obtain the same equation (38) except for different E-terms perhaps ; hence 
the same necessary form for X u , X IK as before. Thus we have the necessary 
form for the entire \-plane. 

The final fact of unique existence is demonstrated as before unless 

f f 

either arg g n = j ^ or arg g in = 

when the values p u , p m may come indefinitely near to the bounding ray 
arg p = (Z + 1 ) ir In . By considering first the case 

pM =pA x ) p«( x ) = °» 

when the solutions y. are known, and by using continuity considerations, the 
fact of unique existence can be established for this case also. 

§ 5. Nature of the expansion. 

In this section we give a notion of the character of w 4 , v. and of the expan- 
sion (7). 

Let us first develop the character of u u in the case n = 2p, — l. Clearly 

(39) w K =c 1 y 1 + ■•• + c n y n 

where we have for the determination of c x , c 2 , • • •, c n the equations 

(40) c l W i (y l )+c 2 W i (y 2 ) + -.. + c n W i (y n ) = (i=l, 2, ■••, »). 

Substitute in the matrix 11^(^)11 the expressions given in (31), (32), (33) 
and remove the factor />*' from the ith column ( i = 1 , 2, • • -, n); write for p 
the values (36). This matrix then takes the form 
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Pn 


P a ■ 


• Ply. « W M -Pi, ,1+1 


... 


e «y P ln 


(41) 


Pn 


P 

J. 22 


■ P e *V P 

■L 2n *> " J- 2, n + 1 








P nl 


P 


I' >+i 






where the P.. are 


of the form 






V 

(42) 




5 ^ 






and 














I 


= 2^v/-T^ 2,rl// - 


Ifti 





The c 4 are proportional to the minors of any row, say the first. The minors 
of the elements 1 to p contain all the factors 

n W{ 

n « ^. 

«=,i+i 

besides a factor like (42), while the minors of the element i, i > fi, contain all 
these except e'' Wi/M > besides a factor like (42). We conclude that 

c>=Qi (i = l,3, ■■-,?), 

(43) ,3 

c i = e *»Q 4 «=0 + l, ? + »,■■•,»), 

where Q ; is of the form (42). 

If we substitute into (21) the values (36), we obtain 



(44) 

where Y is of the form 

(45) 



Vi 



vi- x lx—a\ 

. = e ™» {*-«' Y. 



(» = 1,2, •■•, n), 



&(«) 



^ 

Z" 1 ' 



Placing the expressions for c. and y. of (43) and (44) in (39), we derive the 
form of u u ; likewise the form of u m , v u , v ul is obtained. We have finally 

»=1 i=n+l 

H— 1 „ Wj / 3-6 \ ,„ I x-a\ n Wj I x-a\ 

«ni = X e -^ »- ' CT m + e' ^ »-> f7" II(i + £ e -^ »- ' U m , 

»=1 i=n+l 

»=1 i=,i+l 



(46) 



<=m+i 
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v = Mr v-i + ^Ek- , r = - 2 fcr v — i _ 2 ^-i.?». . 

Here £7, F" are of the form (45) and ra = 2/x — 1 . 
When m = 2/i we get in the same manner 

ftzj ,„. w <('•_-< l, ) i"( x ~") -i"(*-") -^ ;"-"''■ f*~-) 

«ih=Z« w » y "- a 'U lH +e Kl >- a 'U nil -t-e v "- a; ^11,^1+ S e "" Kb ~ a/ U Ui , 

M -l _ Wilx-b\ i>( x ~ a \ V ( x ~"\ n _r "'' ( x ~ a \ 

(47) ^ »„=.£« '*«.^V K +« '^-<" , F I(1 - r -e U - J F I)(l+I + ^e '»^-« ; F It! 

i=M+2 

/ = 2«r V —\ -\ , 6 = 2«r v-H • 

w n 

In both (46) and (47) it is to be noted that the 2-terms are not important 
except at x = a and x = 6 , since the real parts of the exponential terms are 
large and negative for I large and a < x < b . We omit the fuller discussion, 
to be made by (46) and (47) . It is worthy of remark that the simplest case 

n=l, L(u) = u, W l (u) = u(a) — u(b) 

yields the Fourier series. The general expansion is clearly of a similar nature. 

§ 6. Convergence of expansion tof(x). 
We have expressed the sum of n terms of the expansion (7) in the form 

/=— L= f fG(x,s; X)f(s)dsd\, 

ZlT V 1 JtJa. 

where T is some contour which contains X, , X 2 , • • • , X n but no other character- 
istic values within it ; G is explicitly defined [see (10)] in terms of the solu- 
tions y i of L ( u) -f Xu = . By means of the known asymptotic character of 
the y t we determine the character of G , and then evaluate the contour integral 
as r enlarges without bound. 

For convenience let us use the notation [10] for an expression 

a E 

where w, a, • • • are independent of p but need not be constants. 

Trans. Am. Math. Soc. 86 
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Theorem. Let f(x) be made up of a finite number of pieces in the 
interval a = x = b, each real, continuous, and with a continuous derivative. 
Unless x = a or x =5, the expansion for f(x) connected with the differ- 
ential equation L(u) + Xw = and the regular boundary conditions 
W i (u) = 0(i = l,2,.-.,n) converges to [f(x + 0) +f(x - 0)]/2. At 
x = a the series converges to 

«,/(« + °) + « 2 /(&-0) 

and at x = b to 

b i f{a + 0) + b 2 f(b-0), 

where a x , a 2 , b x> b 2 are constants independent off(x).* 

Proof We restrict ourselves first to tha case n = 2//. — 1 . At the end we 
outline a similar proof for n = 2/i . 

Let the contour T be taken as a circle | X I = k in the X-plane. If the trans- 
formation X = p n be made, we find 

/= f C np n - l G(x, s; \)f(s)dsdp + f f np n - 1 Q (x , s ; X )/( s ) ds dp . 

In this expression y x is the segment of a circle \p\ = k' lying on a region S, 

s,-. h^ SvgpS (A±Dz % 

1 n 6r n ' 

and 7 2 is the segment of the same circle lying on the adjacent region S, 

z n ° r n 

We confine our attention to the partial integral 

(48) I"-* = f f np"~ l G(x,s;X)f(s)dsdp, 

taking a<cc<5. It is found to tend toward f(x — 0)/4. By considera- 
tions of symmetry the limiting values of the remaining partial integrals are 
determined. The first part of the theorem then follows. 

Let B be the bisection point of <y x and A , O its end points. Take AB as that 
segment of y : for which 

(49) K(pwj<0 

and B C as that segment of y : for which 

(50) #(/*>„) >0. 

* For definitions of regular conditions see §4. The restrictions imposed on/(i) might be 
lightened but only at the cost of brevity. 
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The indices 1 to n of w i are chosen as in (29) and (30). We must break up 
I a ~" in two parts I a j^ and I^ x , both of which are attacked in the same way. 



To treat G on A B, write 



(51) 



G(x,s; ^) = ±|£i/ i (x)y i (s) 



+ if x > s , 
— if x < s, 



where the y i of (11) have been taken as the y { of (21). Then y t is the coeffi- 
cient of t. in 



j£-»>(,) 2/;- 2 >( S ) ... y^- 2 \s) 



Vii s ) y 2 ( s ) 



y n ( s ) 



<- 2) (s) tf-^e) ••• 2/r 2) («) 
^( s ) 2/ 2 ( s ) ' ••• y„(«) 

But by (21) 

(52) y[P(s) = e <> w >< 3 - a) p J [ w{ ] . 

If these expressions be substituted above we find 

(53) y i ( « ) = «-p-^-> -^ [ - «,, J * 
since — wjn is the coefficient of t. in the expression 




* The denominator contains the factors e pWiC3 a) which we divide into the coriesponding i-th 
columns of the numerator. It is to be remembered that, if t 4= 0, 



[ 



t] [t\- 
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In fact this coefficient m. satisfies the conditions 

I 

n n 

J^ t w{m i =0 (j = 0, 1, ■■-,n — 2) ^w n r l m i = l. 

These show that m. has the stated value. 
By (51) 

n 

G(x,s; ^) = + iT,y i (x)y.(s) 
since a < s <_x for the partial integral 1°^ . Also 

w^G) = -hiw ia (vM( s ) + * £ w»<y,)y,(*)- 

From these equations it is clear that if we multiply the columns 1 to /u. of 

iV(a:, s; X) in (10) by + ^y t (s) (i = 1, 2, • • -, /*), the columns /* + 1 to w 
by — Jy f (s)(i = /u. + l,/ti+2, • • • , w ) and add them to the last column, this 
last column has the elements 

+ I>«("0y<(*)> t.w^yMis)- ± w ia { yj )U s ) c=i,8,-,»). 

By (21), (28), and (53) these are of the forms 

n P i=i 
(54) 

1 / M " \ 

~ s=l( Seo^'-VE— |8rf' +1 ]+ I] e<' w * a - s) pP i [a i w ki+l ] ) (*=1, 2, •••, «). 

In the new notation equations (21), (31), (32), (33), and (34) give also 

y, = «"*■-> [1] (*=1,8,...,»), 

>*<[a.wf] 0=1,2,--..,a«-1), 

/'( [« 4 <] + e^M*-") [£«£] ) (i=^), 

L /<e«<^ [&«*«] (i=K + l, * + &, ■■•, »), 



(55) Fi( y> ) = . 



A (^)=-n/» 4i - II e pw ^- a) ([^ ]+e' >M ^- o) [^]). 

We introduce these expressions into (10) where the elements of the last column 
of JV have been modified to the form (54). The factors of the denominator A ( X ) 
we can divide into the numerator as follows : the factor — 1 into the last column, 
the factor p h into the ( i + l)-th row, the factor e'""' < * _a) into thej-th column, the 
remaining factor A into the /ti-th column. We thus obtain 
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(56) G(x,s; \)=— i=i x 

e pw plx—a) r J "] u , 

r a w *n . . . M V'([ a i]+ e ' wv " [gy J ... f 8 w^l Y e^fi-^B w* I+l l : 

L a i w iJ rfl j +e p^(»-a)r^ j L^i w » J ^ e L^i w i J i 

n 

_ £ e p«'/«-»)[a n 7D*-+ 1 ] 

j=n+i " J i 

From this form of G we can easily determine the limit of I a ^ . 
Assume the sequence of circles so chosen that their distance from the nearest 
point p u , p IU is at least d~> 0.* Then on AB when \p\ is large 

(57) [0 O ] + «*»*»-«> [0J > Jf> 0. 
In fact if 

[«o]+« P "" ( *-" ) [«i]-*. 

where 1 1 is small, we infer at once that the point p lies near to p u or /> m . 

In view of (57), every element of G given by (56) is bounded on AB since 
the exponents have negative or zero real parts [see (29) and (49)]. Further, 
when this expression for G is used in I ^ the factor np n ~ l cancels. 

Consider those terms of I a Xs which do not contain the first element of the 
last column ; s is confined to this column. When we integrate these terms as 
to s, they have the form Mjp (M bounded), as an integration by parts between 
the points of discontinuity of y( s ) shows. But these terms all contain as factor 
an element of the first row, not the last element, beside other bounded elements ; 
since these first row terms are small for | p | large, except the /a-th one when p is near 
to B,-\ we conclude that the p-integration of them will yield only small terms. J 

*This refers to the minimum distance of pu , pin from any point of the circle. In the proof 
of the theorem of \ 4, pu, pm (corresponding to %u, %iu) were found to he approximately equally 
spaced along the bisectors of Si and S 2 . The above construction of a sequence of circles is there- 
fore possible. 

t Then pu>u. is a pure imaginary. 

J'B £ 
Ap dP 



< maximum of e X x- • 
2» 



This is small everywhere on A B if c is small. Even if e is small except for a little part of AB , 
this integral is small. 'Small' means of course indefinitely small as \p\ becomes indefinitely 
large. 
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Accordingly we need only consider the terms which contain the first element 
of the last column. This element has a coefficient 1, as it had originally in 
G(x, s; A.). Here again we can restrict ourselves to terms in the brackets 
which do not contain 1/p. Finally therefore 

I"Ib = I 7— f f f it, e»™<<*-*w)f( s)dsdp + e, 

where e signifies a quantity which tends to as | p | increases. 
Integrating by parts we find 






" f( X — ) € 

e pwix-.) w f( s \d s = -L± 1 + _ (* = 1,2, ■■ ■,/!), 

,-/ K ' p p 



where e is certainly small unless i = p and p is near to B . From this we con- 
clude that 

X B r x 7ri/-^i 
I e" wix ~ s) w.f(s) ds dp = — g , /( x - ) + e , 

and therefore 

(58) I a Is=£ l f(*-0) + e. 

To treat 1"^ we multiply the columns 1 to p — 1 of N{x, s; X) by 
+ i2/t( s )' an( ^ * ne columns p to n by — \y { {s) and add them to the last 
column. The elements of the p-tia. column can be written in this case 

" V L * J '-^^ (t = l,2, --.n), 



e( ,„ v( a-« [ O ] + [ 0, ] ' e*v«-»> [ 6> ] -f [ 0, ] 

and all the elements are bounded as before. The important term again comes 
from the first element of the last column, and we find precisely as before that 

(59) I^ = ^f(x-0) + e. 

Adding (58) and (59) we find 

and by symmetry 

i a y r = if( x - ° ) + e ' J v = if( x + ° ) + e ' 7 ;:* = */(* + ° ) + e - 

Hence we find 

T „. h f(x-0)+f(x+0) 
J f " 2 + 

This proves the first part of the theorem. 

When x = b the preceding work must be modified, since now the elements 
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(i + 1 to n of the first row in (56) have the form [1]. Other elements of 
the last column now become important beside the first. An integration by parts 
with respect to s followed by a ^-integration gives us terms 

b 1 f{a + 0) + b 2 f(b-0). 

The same result must of course hold at x = a. 

When n = 2/j. the attack is almost the same. One obtains a determinant 
expression for G by modifying the last column and distributing the two factors 
of A into the £i-th and (^t + l)-th column.* The elements are bounded on y x 
in this case and one obtains first I"~ r ' which tends toward \f{x — 0). The 
remainder of the proof is as before. 

*See (38) in which A is factorable into 



